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00 . 

<^ , We present a lattice calculation of Liq, one of the low energy constants in Chiral Perturbation 



Theory, and the charged-neutral pion squared mass splitting, using dynamical overlap fermion. 
5_i ■ Exact chiral symmetry of the overlap fermion allows us to reliably extract these quantities from the 



difference of the vacuum polarization functions for vector and axial-vector currents. In the context 
of the technicolor models, these two quantities are read as the ^-parameter and the pseudo-Nambu- 
Goldstone boson mass respectively, and play an important role in discriminating the models from 
others. This calculation can serve as a feasibility study of the lattice techniques for more general 
technicolor gauge theories. 
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Spontaneous chiral symmetry breaking (S^SB) of strongly interacting gauge theory may 
provide a natural mechanism for the electroweak symmetry breaking. A class of new physics 
models based on this idea, so-called the technicolor models, has been studied extensively [l| • 
In most of those models, massless techni-quarks with weak charge are introduced; the weak 
gauge bosons acquire masses from their S^SB. The S'-parameter may then be sizably af- 
fected, for which those models can be strongly constrained through the electroweak precision 
measurements ^]. Another characteristic signal of the technicolor models, that may be ob- 
served at the LHC experiments, is the presence of extra Nambu-Goldstone bosons (NGBs) 
which are not eaten by the weak gauge bosons. They are called the pseudo-NGBs (pNGBs), 
since they must be made massive by introducing explicit breaking of the chiral symmetry of 
the techni-quarks in a model dependent way, otherwise they would remain massless. Since 
the S-parameter and the pNGB mass are consequences of strong dynamics of the underly- 
ing theory, non-perturbative framework is required for their calculation. In previous studies, 
some model was involved in the calculation, e.g. 3]. 

In this work we consider two-flavor QCD as a testing ground of our method and demon- 
strate that the first principles calculation of those quantities are possible. In this context, the 
S'-parameter corresponds to L r 1Q (or l§ in another convention), one of the low-energy constants 
of the chiral perturbation theory (ChPT), as S=-l§it[L r w {u) - {ln(fi 2 /m 2 H ) - l/6}/192vr 2 ] 
with a renormalization scale jx and the Higgs mass m# 2J. L[ is related to a differ- 
ence of vacuum polarization functions between vector and axial- vector currents Hy]^ A (q 2 ) = 
Uy\q 2 ) — Il^(g 2 ) near the zero momentum insertion. (A formula will be given in (JHJ).) 

For the pNGB mass, a mass formula that is valid for a wide range of technicolor models 
and breaking patterns is known 4|. The formula contains a nonperturbative part written 
in terms of the vacuum polarization functions. The charged pions in two-flavor QCD is an 
example of pNGB, as the electromagnetic interaction explicitly breaks SU(2) chiral sym- 
metry and gives a finite mass even in the massless limit of up and down quarks The 
corresponding mass formula is known as the D as- Guralnik- Mat hur- Low- Young (DGMLY) 
sum rule |(J 

3a r° 2 <z 2 ng^fa 2 )U,-o m 
m « ± = -^l dq P ' (1) 

which gives the mass of charged pions at the leading order of the electromagnetic interaction. 
Here / denotes the pion decay constant in the chiral limit. Note that neutral pion is massless 
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in this limit. 

In the continuum theory chiral symmetry guarantees that the difference Uy^_ A (q 2 ) exactly 
vanishes in the absence of both explicit and spontaneous chiral symmetry breaking. Any 
remaining difference in the absence of explicit breaking thus signals the S^SB. Therefore, 
the use of exactly chiral fermion formulation is mandatory in the lattice calculation, in order 
to avoid fake contributions to Uy\ A (q 2 ) due to non-chiral lattice fermion formulations such 
as the Wilson- type fermions. Here we use the overlap fermion j3] , which respects exact chiral 
symmetry at finite lattice spacings. Employing this fermion, we have successfully done a 
precise calculation of the chiral condensate 8| , which also requires excellent chiral symmetry 
to control systematic errors. 

We perform a two-flavor QCD calculation on a 16 3 x 32 lattice at a lattice spacing a = 
0.118(2) fm determined with the Sommer scale r =0.49 fm as an input The quark mass 
in the lattice unit is rh q =am q = 0.015, 0.025, 0.035, and 0.050, which roughly cover the range 
between 1/6 to 1/2 of the strange quark mass. The global topological charge Q is fixed to 
ensure the exact chiral symmetry [9(. The main simulations are done in the Q=Q sector, 
using 10,000 trajectories. At rh q = 0.050, the simulations are also performed in other two 
sectors (Q = —2 and —4) to estimate the finite volume effect due to fixing Q 10] . For each 
sea quark mass, the measurements are made at every 50 trajectories. Statistical errors are 
estimated from a jackknife analysis with 100 jackknife bins each containing two consecutive 
measurements. Details of our configuration generation and the pion spectrum and decay 
constant analysis are found in 9j] and HI], respectively. 

We calculate the current-current correlators for vector and axial-vector currents to 
obtain the corresponding vacuum polarization functions. We use as the vector current 
V/2 —Zqij^l — aD/2m )q2, where q± and q2 represent different flavors of quarks, D the 
overlap-Dirac operator in the massless limit, and m =1.6. The axial- vector current is 
the same but 7 M is replaced by 7^75. The factor (1— a-D/2m ) is necessary to make the V and 
A form an exact multiplet under the axial transformation. Because of this exact symmetry, 
then leading to the strong correlation, even the lattice artifacts and statistical fluctuations 
cancel between VV and AA correlators except for the effects of S^SB. Indeed, the statistical 



errors are much smaller than those of the previous calculations of the VV correlator 12]. 
The common renormalization constant Z = 1.3842(3) is determined nonperturbatively 11]. 
Since the continuous rotational symmetry is violated on the lattice at 0(a 2 ) and the 
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currents we use are not conserved (cf. 13j), the general form of the current-current correlator 
reads 

n^(g) = ^^(oir^c^jrHo)] |o> 

X 

oo oo 

= E B[ j\%f n ^+ E c { r\%? n - l %) 2m -\ (2) 

n=0 ra,m=l 

where J = V or A. Bj and Cj 2 '^ are scalar functions of lattice momentum q^lixn^j L 
with n M an integer ranging from —L/2 + 1 to L/2 (L=1Q or 32 for spatial or temporal direc- 
tion, respectively). In the continuum limit, only and Cj 1,1 ^ survive. B^ could contain 
a power divergent contribution due to a contact term, but the exact symmetry present be- 
tween the vector and axial-vector currents guarantees that this contribution cancels in the 
difference Tlvuv ~ Coefficients other than Bf^ and Cf^ represent lattice artifacts. 

In the difference ITy^ — II^j,, these lattice artifacts are negligible as numerically confirmed 
below. 

We define a measure of the Lorentz-violating lattice artifacts by 

which contains all of B^ and cj™'" 1 " 1 but Bf^ nor Cj . Figure [T] shows Aj for J = V 
and A (top) and their difference (bottom) as a function of q 2 at m g =0.015. While we 
observe statistically significant non-zero values of Aj depending on q 2 , the difference is 
orders of magnitude smaller than the individual Aj. Similar plot is obtained for m g =0.050. 
This indicates that the Lorentz-violating lattice artifacts indeed cancel in the difference 
Tivav ~ ^Auv an d are insensitive to S^SB or m q . Neglecting the Lorentz-violating terms, we 
analyze the difference 

LV^ - U Aflu = {q 2 5^ v - q^ L q v ) 11^ - qpq v U$L At ( 4 ) 

where IIy^ A and Ily^ represent the transverse and longitudinal contributions, respectively. 



First we calculate L r w {ji) from IIy_ A . At the next-to-leading order, ChPT predicts [l4j 



f 2 In + i - H(x) 

41^ 2 ) = J f 2 -slim , (5) 



H(x) = (l+x 
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y/l + x In ^ 



Vl+x+l 



(6) 
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FIG. 1: q 2 dependence of Aj (J = V or A) (top) and their difference (bottom). The result for 



m 9 =0.015 is shown. 



where x = 4m 2 / q 2 , and \i is a renormalization scale set equal to the physical p meson mass 
m p . Using the measured values of fn, and /„-, we fit the data of q 2 Tly\ A at four quark masses 
with (j3j) to obtain L r 10 (m p ) varying fit range of q 2 . Correlation among the data points are 
ignored since each of all the data comes from different sea quark ensemble (also see below). 
It turns out that the fit including only the smallest q 2 point (g 2 =0.038, which corresponds 
to (320 MeV) 2 ) gives an acceptable x 2 /dof (~0.5). The fit is shown in Fig. [2] as a function 
of m q (circles and solid curve). Once the second smallest q 2 (~ (650 MeV) 2 in the physical 
unit) is included the fit becomes unacceptable (x 2 /dof ~ O(40)). This may indicate the 
breakdown of the chiral expansion at such a large q 2 . Our result from the smallest q 2 data 
is L r l0 (m p ) = —5.22(17) x 10 -3 . Here, the error is statistical only. 

We estimate the systematic error due to higher order effects of the chiral expansion using 
a modified fit function to cover a wider range of q 2 (see below) . We obtain a slight negative 
shift, 0.3 x 10 -3 , which is added to the systematic error. The finite size effect may be sizable 
in the pion-loop effects, which is the third term in dSJ), since the lattice volume (1.9 fm) 3 
is not large enough. We estimate its magnitude by replacing the momentum integral with 
a sum. and m n are also corrected following 15(. Taking these corrections into account, 
we fit the data at the smallest q 2 to (jSJ) and obtain L\ Q (m p )\v=oo — —5.74(17) x 10 -3 with 
X 2 /dof= 2.3 as shown in Fig. [2] (triangles and dashed curve). We take the difference between 
these two results as an estimate of the systematic errors. We then quote 



-5.2(2)(±§)(±g) x 10 



(7) 



where the first error is statistical, and the second and third are the estimated systematic 
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FIG. 2: q 2 n ( y ] _ 

A I q 2 =0.038 as a function of quark mass. The fit results with (solid) and without 
(dashed) finite volume correction are shown. 

uncertainties due to higher order effects in q 2 and the finite size effect, respectively. Since 
only one value of q 2 is included in the fit, the error from the chiral fit may be underestimated. 
Furthermore, other sources of uncertainty, e.g. finite lattice spacing and lack of a dynamical 
strange quark, exist. Nevertheless, ([7]) is already consistent with the experimental value 
-5.09(47) x 10" 3 Q. 

Next, we consider the squared- mass splitting between charged and neutral pions. The 
splitting in the chiral limit solely comes through the electromagnetic interaction and is writ- 
ten by the integral of q 2 Tly\ A as given in ([T|). In order to avoid possibly large discretization 
effects in the large q 2 region, we separate the whole integral region into two parts at g 2 =2.0, 
and estimate each part as follows. 

For the lower q 2 region (< 2.0), we fit the data to an ansatz 



fU^ A \q 2 ) = -fl 



q 2 f 2 v q 2 fl <f 



X{q 2 ) 



q 2 + rhy q 2 + m\ 24tt 2 1 + x 5 (Q 



2\4 ' 



(8) 



where Q 2 = q 2 /m 2 with rh p the physical p meson mass in lattice unit. Here and in the 
following Xi denotes a fit parameter. We introduce poles of the lowest-lying state for both 
vector and axial- vector channels with masses my^ and decay constants fv,A- We put the 
constraints f 2 = fv~ f\ and f ' attla = fyrhy among them so that they satisfy the first and 



second Weinberg sum rules 
and my = X2 + 24 m 2 . The function X(q 2 ) is either 

ml\ 1 



171 ]. We also assume a linear dependence on rh 2 : fy = X1+X3 m 2 . 



In 



mi 



-H(4m 2 Jq 2 )+x 6 Q 



or x 6 Q 2 MQD- 



(9) 
(10) 




FIG. 3: The fit results with ([9|) (dashed curves) and (jlOp (solid curves). The results in the chiral 
limit are also shown. The statistical errors shown are smaller than the symbol size. 



Then, the function (jSJ) behaves as 0(q 6 ,q 6 In q 2 ) at large q 2 in the chiral limit, which is con- 
sistent with the asymptotic scaling predicted by the operator product expansion (OPE) 18]. 
Taking © for X(q 2 ), IT^^g 2 ) reduces to the ChPT prediction © when Q 2 p < 1, while 
(TTU|) gives a logarithmic term in the large Q 2 region as expected by OPE. 

We fit the data at q 2 < 2.0 using the measured values of and rh n as shown in Fig. [31 
We have only attempted an uncorrelated fit since the full covariance matrix is likely ill- 
determined for many data points and free parameters in this fit. Both (Q and (fTUl) fit the 
data quite well and indeed give a reasonable x 2 /dof, though the latter is slightly better. 
Integrating over q 2 in the chiral limit, we obtain m 2 ±|g2< 2 .o= 676(50) and 811(12) MeV 2 for 
and ffTUl) . respectively. The difference arises from the chiral extrapolation around g 2 =0.1- 
0.2, since ([9]) contains the chiral logarithmic term. Recalling that in the determination of 
L r 10 the ChPT formula fits the data only at the smallest q 2 and (fTUl) fits the data better than 
(Q, we take the central value from the fit with ffTUj) and the difference as a systematic error 
due to the chiral extrapolation. 

Expanding g 2 !!^^ around q 2 = in the chiral limit and comparing with (jSJ), we obtain 
L r 10 (m p )= -f 2 (2x 2 - / 2 )/(8x 2 x 2 ). With the fit results for ©, this gives L r 10 {m p ) = -4.9 x 
10 -3 . The difference from the central value is added to the systematic error from the higher 
order effect, and included in (|7j) as already mentioned. 

The remaining part of the integral (q 2 > 2.0) is estimated based on the OPE, which 
^ a 6/((/ 2 ) 3 in the chiral limit for large q 2 up to a logarithmic term. 



predicts Hy_ A (q 



Assuming Ily^l^ =0 = a 6/(<? 2 ) 3 & t 9 2 =2, the fit result with (fTUl) gives a 6 =— 0.0035. In the 



7 



estimate of the final result, we use a phenomenological value in the range [—0.001, —0.01] 



GeV 6 



19] to be conservative. An integral then gives fn 2 r ±\q^>2.o = 182 (149) MeV 2 . 



Summing up the two parts, we obtain 

m 2 ± = 993(r2)(± 13 5 )(149)MeV 2 , (11) 

as the pion squared- mass splitting in the chiral limit. The first error is statistical; the second 
and third ones are due to the chiral extrapolation and the uncertainty in a 6 . The result is 
reasonably consistent with the experimental value at the physical quark mass [1261 MeV 2 ]. 
In addition to the errors quantified above, other sources of systematic errors may still remain. 
We do not expect, however, substantial systematic errors other than those estimated above, 
since the integral is dominated by the q 2 ~ region where the integrand q 2 Uy\ A / f 2 in the 
chiral limit is strongly constrained by the first Weinberg sum rule [q 2 Hy\ A ]q2 =0 / f 2 =l. 

In this letter, we have demonstrated that the S'-parameter and the pNGB mass can be 
calculated using the lattice QCD technique. Since these quantities are generated solely 
through S^SB, the exact chiral symmetry on the lattice plays an essential role to prohibit 
contaminations from the explicit breaking. The method is general and the application to 
other vector-like gauge theories with arbitrary number of colors and flavors is straightfor- 
ward. Thus with this method the lattice technique is able to directly investigate physical 
quantities relevant for the LHC phenomenology. In addition to these quantities, we can also 
calculate a§ and the strong coupling constant using the data in the large q 2 region. The 
results will be reported in a subsequent paper. 
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